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We investigate preheating of an inflaton field (f) coupled 
nonminimally to a spacetime curvature by the method of 
Hartree approximations. In the case of a self-coupling inflaton 
potential V((j)) = A(/>*/4, the dynamics of preheating changes 
by the effect of the negative £^. We find that the nonminimal 
coupling works in two ways. First, since the initial value of in- 
flaton field for reheating becomes smaller with the increase of 
1^1, the evolution of the infiaton quanta is delayed for fixed A. 
Second, the oscillation of the infiaton field is modified and the 
nonadiabatic change around = occurs significantly. That 
makes the resonant band of the fluctuation field wider. Espe- 
cially for strong coupling regimes |^| 3> 1, the growth of the 
inflaton fluctuation is dominated by the resonance due to the 
nonminimal coupling, which leads to the significant enhance- 
ment of low momentum modes. Although the final variance 
of the inflaton fluctuation does not change significantly com- 
pared with the minimally coupled case, we have found that 
the energy transfer from the homogeneous inflaton to created 
particles efficiently occurs for ^ < — 60. 
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I. INTRODUCTION 

Recently it has been recognized that a post- 
inflationary reheating begins by an explosive particle pro- 
duction which is called preheating [0,p[- As compared 
vifith the old theory of reheating Q , preheating is a non- 
perturbative process by which the energy of the coher- 
ently oscillating inflaton field is efficiently transferred 
to scalar or other particles by parametric resonance. The 
existence of this stage would drastically change the ther- 
mal history of the Universe, and many authors have in- 
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vestigated several issues such as baryogenesis at the GUT 
scale [Q, the non-thermal phase transition (5|, the for- 
mation of topological defects fpl , and gravitational waves 

S|. 

There are two typical models of preheating in the 
chaotic inflation scenario. One of which is the massive 
inflaton w?(\? jl case. In this model, another scalar field 
X coupled to the inflaton field can be resonantly ampli- 
fied. There are several numerical works by making use 
of Hartree approximations B and fully nonlinear lattice 
simulations po|jll| ]. At the linear stage of preheating, 
the evolution of produced fluctuations can be analyzed 
analytically by making use of the Mathieu equation [Q . 
When the coupling constant is large, parametric reso- 
nance turns on from the broad resonance regime, and an 
efficient particle production is possible overcoming the di- 
luting effect by the expansion of the Universe. However, 
if the inflaton field (/) does not couple to another scalar 
field Xi the inflaton fluctuation does not grow because 
the resonant terms are absent in a single field model. 

Another model is the self-coupling inflaton A(/()^/4 case. 
In this case, the fluctuation of the inflaton fleld is gen- 
erated even if we do not introduce another scalar field 
coupled to inflaton. The non-equilibrium scalar field dy- 
namics are studied by fully nonlinear calculations p3[ , 
and are also investigated by making use of Hartree or 
large- iV approximations M-ttSJ. Since the Ricci scalar 
averaged over several oscillations of the inflaton field van- 
ishes, the system can be reduced to the theory of the 
Minkowski spacetime. Making use of the properties of 
the Lame equation, we can find that the growth rate 
of the inflaton quanta is not as large as that of the x 
field in the w?(f? 12 -I- g^(tP'x^ /2 model. However, the in- 
flaton fluctuation grows quasi-exponentially and reaches 
the maximum value {5(1)^) ^ lO^^Mp^. In both models, 
the backreaction effect on the oscillating inflaton field 
plays an important role to shut off the resonance. Also, 
a rescattering effect becomes significant at the final stage 
of preheating. 

Generally, it has been assumed that a nonminimal cou- 
pling of the inflaton field or other scalar fields can be neg- 



ligible in preheating phase. However, it was first pointed 
out in the massive inflaton model that the nonminimally 
coupled X field can grow by the oscillating scalar curva- 
ture |19| especially when the coupling constant ^ is neg- 
ative. Later we investigated the structure of resonance 
including the (j)-x coupling where x is nonminimally cou- 
pled to the scalar curvature [Q. In that model, reso- 
nance occurs in the wider range as compared with the 
minimally coupled rn^cjp /2 + g^cfp"}^ /2 model, and the 
final variance of x-p£^rticles can be larger if ^ < — 4. In 
other inflation models such as the higher-curvature grav- 
ity, x-particles can be produced quite efficiently when |^| 
is of order unity Q . 

By these works, we can expect that the nonminimal 
coupling modifies the structure of the resonance also in 
the infiation model without other fields and makes pre- 
heating more efficient. In the massive inflaton model, 
since ^ is constrained to be small as |^| < 10"'^ by the 
observation of the Cosmic Background Explorer (COBE) 
satellite [^ , the evolution of the system will not change 
so much even if the nonminimal coupling is taken into ac- 
count. On the other hand, in the A0^/4 model, the con- 
straint of 5 is absent for the negative coupling. Especially 
when ^ is lar gely negative, the fine-tuning problem of A 
is relaxed ||2^,Q. It was pointed out in Ref. ||lj,|l^ that 
the nonminimal coupling with the inflaton field would not 
affect the evolution of the inflaton fluctuation because the 
Ricci scalar averaged over several oscillations of the in- 
flaton vanishes in reheating phase. However, since the 
Ricci scalar oscillates by the interaction with the homo- 
geneous inflaton field, there is a possibility that the dy- 
namics of preheating would be altered if we consider the 
nonminimal coupling. In the minimally coupled case, it 
is known that the energy transfer from the homogeneous 
infiaton to produced particles is not efficient. In the non- 
minimally coupled case, however, there is a possibility 
that enhanced momentum modes become wider and the 
efficient particle production is realized. It is worth in- 
vestigating the dynamics of the nonminimally coupled 
inflaton field, because this would affect the non-thermal 
phase transition la] , and the formation of topological de- 
fects l^. 

In this paper, we consider the dynamics of the nonmini- 
mally coupled inflaton field in the chaotic inflation model. 
We suppose that there are no fields other than infiaton. 
The Hartree approximation is used to estimate the vari- 
ance of produced fluctuations. We organize our paper as 
follows. In the next section, we present basic equations of 
the background and scalar fields. In Sec. HI, we study the 
dynamics of the nonminimally coupled massless inflaton 
in the inflationary phase, and present the initial values of 
the inflaton when the reheating stage starts. In Sec. IV, 
we investigate how the nonminimal coupling changes the 
evolution of the system in preheating phase. The final 
variance of the fluctuation is also presented. We give our 
conclusion and discussion in the final section. 



II. BASIC EQUATIONS 

We investigate a model where an inflaton field 
nonminimally coupled with a scalar curvature R: 



l^^->^^' 



y(</)) - ie^02 



(2.1) 



where k^ /Sir = G = Mp^ is Newton's gravitational con- 
stant, and ^ is a coupling constant. V{(l)) is a potential 
of the inflaton field which is described by 



y(0) = -to202 _,_ x^4^ 



(2.2) 



In this paper, we mainly consider the self-coupling infla- 
ton model. 

In the minimally coupled case, the value of the infla- 
ton field at the time tj when the inflation turns on is 
required (j){ti) > 5MpL in order to solve the cosmological 
puzzles. Also, the primordial density perturbation ob- 
served by the COBE satellite constrains the coupling of 
the inflaton fleld as m < 10~^AfpL for the massive infla- 
ton model, and A < 10^^^ for the self-coupling inflaton 
model respectively. 

However, the dynamics of the Universe changes if the 
nonminimal coupling is taken into account. Futamase 
and one of the present authors |2^ showed that in the 
case of the positive coupling, ^ is constrained to be small 
as ^ < lO^'' in both chaotic models. This is due to the 
constraint that the effective gravitational constant 



Ge 



G 



1 - cj^ycj^i 



with 



SttC 



(2.3) 



must be positive to connect to our present universe, 
which yields 4> < 4>c — Mpl/VSttJ. The requirement 
that (pc > 4>{ii) leads the constraint: ^ < lO^''. 



The Lagrangian in Eq. (2.1) is reduced to the system 
in the Einstein frame by a conformal transformation 



— q2 



(2.4) 



where $7^ = 1 — (^K^cfP. Then we obtain the following 
equivalent Lagrangian: 
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and 
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2 _ 1 - (1 - em^'^ 
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(1-^^2^2)2 ' 



yy9)- (^ _ ^^202)2- 



(2.5) 



(2.6) 



(2.7) 



Introducing a new scalar fleld tp ^ 



^ ^ / F(0)d</., 



(2.8) 



the Lagrangian in the new frame is reduced to the canon- 
ical form: 



ds' 



-dt^ 



£ = 



^^-i(vvf 



2k 



vw 



(2.9) 



When ^ is negative, the effective potential ( [2.7| ) shows a 
different characteristic from the minimally coupled case. 
In the massive inflaton ca se, sin ce the potential has a 
maximum at 4> = Mpi,/y^8n\^\ (See Fig. 1), ^ is con- 
strained to be small as |^| < 10~^. As long as we do 
not introduce another field coupled to inflaton, the fluc- 
tuation of the inflaton quanta is not generated in the 
minimally coupled case. Hence the situation does not 
change even if we take into account such small values of 

ie|. 

However, such a constraint is absent in the case of the 
A(/)^/4 model, because the effective potential of inflaton 
has a plateau to lead a sufficient inflation (Fig. f). Es- 
pecially when the coupling is strong as |^| > 1, several 
authors p3|-p5[ investigated the dynamics of the nonmin- 
imally coupled inflaton field in the inflationary phase. 
Hereafter, we consider the model of self-coupling poten- 
tial V{(l)) = X(p'^/A with negative f . 

Let us obtain the basic equations in the self-coupling 
model with the nonminimal coupling. We can derive the 



equivalent equations in both frames of (2.1) and 
However, the equations in both frames are different each 
other if we take the mean field approximation in each 
frame. Generally, it is convenient to consider the evo- 
lution of the system in the Einstein frame. However, 
since the conformal factor is included in the denomina- 



tor of the potential (2.7), we can not replace the mean 



field value of (1/(1 - ^K^ry) as 1/(1 - ^K^{r)r. If 
these replacements were done, we need further assump- 
tions of the derivative terms of the (f) field in order to 
coincide with the basic equations of which the mean field 
approximation is taken in the original frame. Although 
we expect that these derivative terms do not affect the 
preheating dynamics significantly, there still remains the 
ambiguity whether these assumptions are always valid. 

Moreover, the conformal factor fl^ includes quantum 
variable (fp since we consider the quantum fluctuation of 
the 4> field. The metric in the Einstein frame is perturbed 
by making a conformal transformation even if the met- 
ric in the original frame is spatially homogeneous. The 
analysis of other preheating models including the metric 
perturbation are performed by several authors ^^, but 
in the present model, there remains the difficulty of how 
we take the mean field approximation in the Einstein 
frame. In order to avoid these problems which are ac- 
companied by a conformal transformation, we consider 
the basic equations in the original frame and take the 
mean field approximation. 

We adopt the flat Friedmann- Robertson- Walker metric 
as the background spacetime; 



a2(i)dx2. 



(2.10) 



where a{t) is the scale factor. Although we do not con- 
sider the metric perturbation in this paper, we should 
include this effect for a complete study of preheating. It 
is under consideration. 

We obtain the following field equations from the La- 
grangian (2.1), 



^ l^ ^ G,,, = (1 - 2e)V^0V,0 - Q - 2^ 5m-(V0)' 

(2.11) 



-9t..V{(j,) + 2^0(3^. □ - V^V, 



^R(j) - y,0 = 0, 



(2.12) 



where D and V,^ are defined as 

n = dii,{\r-gg^^dv)l'J--g-, V„j,= dV/dcj) respectively. 

Since we consider the quantum fluctuation of inflaton, 
the (j) field is represented by the homogeneous and fluc- 
tuational parts as 



</)(i,x)=(^o(i)+50(t,x) 



(2.13) 



As for the fluctuation part, we impose the tadpole con- 
dition 



(50(t,x))=O, 



(2.14) 



where (• • •) denotes the expectation value. We also make 
the Hartree factorization as 

{5<pf ^ 3((<5^)2)(^0), (2.15) 

{54,f ^ 6((<5^)2)(^0)2 - 3((5</>)2)2. (2.16) 

Then we obtain the Hamiltonian constraint equation by 
taking the mean field approximation in Eq. (2.11), 



H^ = 



3(1 -a) 



\{^l + {5^^)) 



^ - 2e ) (50'^^ 



1 



+-\{4>l + Q4>l{ 



'?) 



+2^{3if (</)o0o + {HH)) - {HH")} 



(2.17) 



where a dot and a prime denote the derivative with re- 
spect to time and space coordinates, respectively, and 
a = iK^{4>'^), H = d/a. 

Next, we consider the equations of the homogeneous 
and fiuctuational parts of the inflaton field. Let us first 
obtain the expression of the scalar curvature R. Since 
R can be written by i? = —g^G^,,^, we easily find the 
following relation from Eq. (2.11) as 



l^^K^ 



-R={1- 6O(V0)2 + W{(f>) - 6^<j>n<j). (2.18) 



Then R is expressed as 



R = 






1 — a 

+ {S(P64>) + 3H{(j>o4>o + {5(t>54>)) - {HH")} 



(2.19) 



Taking the mean field approximation of Eq. ( p.l2| ), tlie 
(j)Q and 5(j) fields obey the following equations of motion, 

ij)o + 3H^o + XM<Pl + 3{Sq^^)) + ^R(bo = 0, (2.20) 



64) + 3HS<p - d^d\5(j)) + {3A(0g + ((50^)) + ^R}5(t) ^ 0. 

(2.21) 

Note that both evolutions of the (/)o and 54) fields can be 
altered by the effect of the nonminimal coupling. 

In order to study a quantum particle creation, we ex- 
pand 54> as 



1 



(27r)3/2 



akS4>k{t)e 



ikx , „t SJ,*/'j-\„ikx\ t3 



(2.22) 



-,t 



where Ofe and aj. are the annihilation and creation op- 
erators respectively. Then, we find that S4>k obeys the 
following equation of motion 



h + 3HS(j)k 



- + 3X{4>l + {S4^^))+^R 



k = o. 

(2.23) 



The expectation values of S(f>'^ and 0^ are expressed as |Tj] 

{34)^) ^ ^ J e\Scf)k\'dk, (2.24) 



(2.25) 



Introducing a conformal time rj = J a^^dt and new 
s calar fields ipo = a4>o and S(pk = a,S4)k, Eqs. ( 2.2C| ) and 
( ^.21 ) are rewritten by 



^ + A^o(^^ + 3a^(502))+(ei?a^-i0 1^0 = 



-f^^'Pk +i^lS^k = 0, 



with 



ujl = k^ + 3X{(pl + a^{54>'^)) + e^c 



Id^a 
a drf' 



(2.27) 



(2.28) 



In the case of the minimal co uplin g, the 3A(po term in 
the time dependent frequency (2.28) causes the resonant 



growth of the inflaton fluctuation. In the nonminimally 
coupled case, £_Ra? term also contributes to the para- 
metric resonance. In this paper, we mainly investigate 
how the inflaton fluctuation ((50^) grows by taking into 
account the effect of the nonminimal coupling. 

As for the initial conditions of the fluctuation, we 
choose the state that corresponds to the conformal vac- 
uum as 



<5<^fc(0) = 



1 






(2.29) 
(2.30) 



We investigate the evolution of the inflaton quanta with 
those initial conditions as a semiclassical problem. Be- 
fore proceeding the particle creation in preheating phase, 
we summarize the evolution of the scale factor and the 
inflaton fleld in the inflationary period in the case of the 
strong coupling in the next section. 



III. THE DYNAMICS OF INFLATION WITH 
THE NONMINIMAL COUPLING 

In this section, we analyze how inflation takes place in 
the presence of the nonminimal coupling term. In the 
inflationary stage, since the fluctuation of the inflaton 
field is small compared with the homogeneous part, we 
can approximately write Eqs. (2.17), (2.18), and (2.2C) 
as follows 



H^^ 



3(1 



2^ 



I + -\4>l + 6ei/0o</'o 



(3.1) 



R 



1 — a 



^(l-6e)0^ + A04-f6^0o(0o + 3i/0o) , (3.2) 



io + 3i?0o 

^2 j,2 



\4)l 



where a — £,K,'^4>t- ^y Eq. (|3 
eliminated to give 



3-ff</)o 



e^^(l-6O0^ 

1 - (1 - 6C)a 1 - (1 - 60a 



\4)l 



^R4>a - 0, (3.3) 

, the scalar curvature is 



00 = 0. (3.4) 



Imposing the slow roll approximations |(/'o/0ol *^ H^ 
(2-26) \4)o/4)o\ < H, and \4)l < F(0), Eqs. (|^) and (jsj) 
can be written as 



H^ = 



t^\4)l 

12(1 -a 



1 - 



se 



1 - (1 - 60a 



iH4po 



\4)l 



l-(l-60a" 



(3.5) 



(3.6) 



In order to see the effect of the nonminimal coupling, 
we consider the stron g n egative coupling case a <C — 1. 
Then Eqs. (|3.5[) and (3.6) are simplified as 



H' 



3H^o 



12a ' 



A03 



(3.7) 



(3.8) 



(l-60«' 
Combining Eqs. (3.7) and (^^), we find the following 



relation 



se 



-H. 



1-6^ 
This is easily integrated to give 



a = aj 



1 



8C , fa 



(3.9) 



(3.10) 



where the subscript / denotes the values when the infla- 
tionary period starts. We find by this relation that |a| 
decreases with the expansion of the Universe. Since the 
scale factor satisfies the relation 



—Xaj 



1- 



se 



ln(^ 



(1 - 6^)q:/ \ai 
we find that a is expressed as 

a = aj exp [Hit - ^{Hit)'^] , 
where 

' 2^ 



1/2 



Hr = 



—Xaj 



12^2^ 



2^2' 



7 



(1 - 60a/ 



(3.11) 



(3.12) 



(3.13) 



At the first stage of inflation, since the 7 term is negli- 
gible, the Universe evolves exponentially. With the pas- 
sage of time, however, 7 term becomes important and 
the expansion rate of the Universe slows down. Fakir 
and Unruh |g^ showed that the inflationary period ends 
when |a| drops down to unity. This can be easily seen 
that the term 



Hcj) 


= 


a 
2Ha 


= 



4e 



(l-60a 



(3.14) 



becomes of order unity when \a\ w 1 in the case of ^ < — 
1. 

This criterion for the end of the inflation can also be 
obtained from the effective potential (^J) in the Einstein 
frame. The slow-roll parameter of the flrst order in the 
Einstein frame is defined by the potential V{ip) as 



e = 



1 (v'WX 

2k^ U(V)/ 



Making use of Eqs. (2.C)-(2.8), this can be written by 

'=«202[i_(i_6Oe«2^2]- (3-16) 



Setting e = 1, the value of a when the inflationary period 
ends is estimated by 



\a.F\ = 



V(l- 240(1 -80 
2(1-60 



(3.17) 



where the subscript F denotes the value at the time tp 
when the inflation ends. Note that JofI = 1 for ^ = — 1. 
jcki?! slowly increases with the decrease of ^ for ^ < 0. 
However, since the limiting value for ^ -^ —00 is |q;f| — 
1.15, inflation ends at \aF\ ~ 1 for ^ < — 1. 

As for the flnal values of (t>Q for the end of inflation, we 



make use of the ones which are determined by Eq. (3.17) 
as 



4>o{tF) = 



v/(l- 240(1 -80 -1 
1671(1-6010 



-,1/2 



MpL. (3.18) 



In the case of the minimal couphng, (poitF) — 0.56Mpl. 
With the increase of |0, <j>o{tF) decreases. We show the 
values of \aF\ and (jioi^F) in various cases in TABLE I. 
Since the large |0 makes the inflaton potential flatter, 
inflation takes place for rather small values of 0o a-s com- 
pared with the minimal coupling case. As for the time 
ti at the start of reheating, we assume ti = tp- Hence, 
the initial value of the inflaton flel d for the preheating is 
^o(ii) = 4'o{iF) estimated by Eq. (B.IJ)- In the next sec- 
tion we investigate the preheating dynamics of the non- 
minimally coupled inflaton field. 



IV. PREHEATING OF THE NONMINIMALLY 
COUPLED INFLATON FIELD 

In this section, we briefly review the minimally coupled 
inflaton case , and after that, investigate the case that the 
nonminimal coupling is included. 



i. ^ = case 

In the minimally coupled case, the evolution of the 
inflaton quanta can be obtained analytically by using 
the stability and instability chart of the Lame equa- 
tion. First, making use of the time-averaged relation 
\{(f^)T — 2(F((/)))t, the evolution of the scale factor is 
approximately written as a ~ 77. Since the scalar curva- 
ture is expressed as i? = ^ ^-| , the time-averaged value 
of R vanishes in the case of the minimal coupling. 

Introducing a dimensionless time variable x = 
\/AA/pl77 and scalar fields / = yp/M pL, ((5^^) = 
(3-15) {5(1)'^) /Ml^, {5(p'^) = a'^{6^'^}, Eq. 0^ is reduced to 
the equation for the theory A0^/4 in Minkowski space- 
time: 



dx'^ 



.{S^') 



(4.1) 



At the stage when the back reaction effect is neghgible, 
the solution of Eq. ( 4J ) can be handled analytically. In- 
tegrating the equation (ff/dx^ + f'^ = 0, we find 



if - fL) 



(4.2) 



where /co is the value when the tpo field begins to oscillate 
coherently. The solution of this equation is expressed by 
the elliptic cosine, 



/ JcO en I JcO\^ ^Coj; /— 

This function is well approximated as 



/ « /co cos(0.8472/eox) 
= /coCos(0.8472%/A(^n(ico)?7), 



(4.3) 



(4.4) 



where we dropped the initial value of x^o for simplicity. 
Note that the frequency of the oscillation depends on 
the values of A and (^o(^co)- As we have showed, initial 
values of the inflaton for preheating changes if we take 
into account the effect of the nonminimal coupling. We 
discuss this issue in the next section. 

Since the inflaton field oscillates periodically, the in- 
flato n qua nta can be produced by parametric resonance. 
Eq. ( 2.27 ) can be written as 



dx"^ 



5fk 



ip\l 



P 



Sifik = 0, (4.5) 



whe re p ^ = k"^ / (XMp-^) is the normalized momentum. 
Eq. ( |4.5| ) belongs to the class of Lame equation when the 
inflaton quanta are not produced significantly. The struc- 
tures of resonance are studied precisely in Ref . ||l6|,|l8| . 
Making use of the chart of stability and instability bands 
of the Lame equation, the growing mode of Sipk exists in 
the narrow range: 



lp<p'<V3f^ 



(4.6) 



where variables with a tilde ^ means the amplitude of the 
variables in this paper. 

As compared with the model of the massive infiaton 
plus another scalar field x^ "^^ found that the modes 
whose momenta are close to fc = are not enhanced 
and that the maximum growth rate of Sipk is not so 
large. In spite of the small growth rate, {6(fP) grows 
quasi-exponentially by parametric resonance. We show 
in Fig. 2 the evolution of {S(p'^) for the initial value of 
4'oi'ti) = 0.56AfpL, and the coupling A = 10~^^ which is 
constrained by the COBE data [Em. At the initial stage, 
((5(/)^) decreases because the expansion rate of the Uni- 
verse surpasses the particle creation rate. Our numerical 
calculation based on the mean field approximation re- 
veals that (Scfy^) begins to increase after x > 100, and 
reaches the maximum value. 



)f 



1.5 X 10^ 



JWpL. 



(4.7) 



The subscript / denotes the value at the time tf when 
the (50^) reaches the maximum value. This result agrees 
well with the fully non-linear calculation performed by 
Khlebnikov and Tkachev |]l3[. Although our mean field 
approximation does not include the rescattering effect 
(i.e. the mode-mode coupling of the inflaton fluctuation), 
the flnal variance is consistent with the lattice simula- 
tion. This implies that the back reaction effect due to the 
growth of the inflaton quanta plays an important role for 
the termination of parametric resonance. As the inflaton 
fluctuation grows, the back reaction effect prevents the 
coherent oscillation of the ipQ fiel d as is found by Eq. (4.1). 



The resonance band in Eq. (4.5) is shifted away from the 
original position as particles are produced. According to 
the analytic estimation by Greene et al. |jl^, the growth 
of the inflaton quanta stops when {Scf)^) grows up to 






O.O50g. 



(4.8) 



Our numerical value of the termination point is about 
440 and (j)o{xf) « 1.4 x 10~^_Mpl, which yields 



Xf 



)f 



IO-^A/Il by Eq. (|4 



This agrees 
The fl- 



well with the numerical value of Eq. (4.7) 
nal fluctuation in the self-coupling inflation model is 
rather small compared with other preheating models as 
to202/2 + g^cjj^x^f^ or m^(t)^/2 + £,Rx^ /2 ^^. In the 
massive infiaton model, preheating starts from the broad 
resonance regime in certain ranges of the coupling con- 
stant. Hence in this case, the growth rate of the fiuc- 
tuation can take larger values, and the modes which are 
close to fc = are mainly enhanced. In the present model, 
however, the instability bands are restricted as Eq. ([l.(^), 
and parametric resonance is not so efficient. Parametric 
resonance ends when the dispersion of the fluctuation be- 
comes about 20% of the amplitude of the inflaton field. 
We mainly concern that how the effect of the nonminimal 
coupling would change the growth of the fiuctuation and 
the final variance. In the next subsection, we investigate 
this issue precisely. 



^. 5 / case 

Next we proceed to the nonminimally coupled case. 
Its effect can be seen in several places and changes the 
dynamics of preheating. As we have already mentioned, 
the nonminimal coupling alters the initial value of the 
infiaton field 0o when the reheating period starts. With 
the increase of |^|, this initial value decreases and the fre- 
quency of the inflaton oscillation becomes smaller. Sec- 



ond, as is found by Eqs. (2.2C) and (2.21), the nonmin- 
imal coupling affects the evolutions of the ^o and 5(j)k 
flelds through the oscillation of the scalar curvature. At 
the stage before the inflaton quanta do not grow signifi- 
cantly, the scalar curvature is approximately written by 



R. 



1 - (1 - 60a 



(4.9) 



where we used Eqs. ( |3.2| ) and (3.3). Note that the oscil- 
lation of the scalar curvature depends on ^ significantly. 
Hereafter, we investigate how the nonminimal coupling 
would have influence on the evolution of 0o a-nd 5(1) k fields. 



We can rewrite Eqs. (2.26) and (2.27) by taking notice 
of the relation R 



6 d^ a 
o^ diq'-^ 






/M 1 + 3 



dx^ 



5fk 



{Sv' 



p2 +3/^1 



a'^Rf = 0, (4.10) 



^-q]^'R 



P 

Stpk = 0, 



(4.11) 



where R = R/{XAIpjJ. In the conformal coupling case 
(^ = 1/6), although it breaks the constraint ^ < 10~^, the 
system is reduced exactly to the equations in Minkowski 
spacetime. 

In the case of t he m inimal coupling, the last terms 
in Eqs. (4.10) and (4.11) become negligible after several 
oscillations of the (po field. In the positive ^ case, ^ is 
constrained to be ^ < 10"'^ for the infiationary period 
to proceed well. We have numerically checked that the 
evolution of the system is similar to that in the minimally 
coupled case. Hereafter, we concentrate on the case when 
^ is negative and fix the self coupling constant as A = 
10~^^. At the end of this section, we will briefly comment 
the case when A is not fixed. 

In the case of— 0.1 < £, < 0, since the initial value of 
I a I is much smaller than unity, the effect of the nonmin- 
imal coupling does not alter the evolution of the infiaton 
significantly compared with the minimally coupled case. 
When —10 < £, ^ — 0.1, ^ effects begin to be important 
and change the evolution of the system. In this case, it 
appears mainly in the value of (j^oiU) and 0o(ico)- Let 
us consider the case of ^ = —1 as an example. In this 
case, preheating begins when (j)o{ti) — 0.1995Mpl and 
\ai\ — 1. Since the value of (/'o(^co) is small compared 
with the minimally coupled case, the period of the os- 
cillation of the inflaton field becomes about five times 
longer than the ^ = case (See Fig. 3). The growth of 
the inflaton quanta is delayed due to this effect. Since |a| 
is of order unity at the end of inflation, we expect that 
the evolution of the inflaton field and the inflaton quanta 
would b e affe cted b y the presence of the (^— h)a^R terms 
in Eqs. (4.10) and (4.11). However, after several oscilla- 
tions of the inflaton field, (po drops down under 0.02AfpL, 
which means that |a| rapidly becomes much smaller than 
unity. We have numerically checked that these terms 



and 



are negligible compared with the /'^ I 1 + 3 ^-^' 

3/^ ( 1 + S ) terms for x > 100. Hence the structure 



of resonance for x > 100 is almost the same as that in 
the minimally coupled case. 

Since the particle creation rate is smaller than the ex- 
pansion rate of the Universe at the initial stage of pre- 
heating, (50^) decreases at the first stage (See Fig. 2). 
The numerical calculation shows that ((50^) begins to 
grow after x w 400, which means that the evolution 
of the inflaton quanta is not affected by the £, effect 
directly. Although the growth of {d(j)'^) is delayed due 
to the decrease of the initial value of the inflaton field, 
the final variance is almost the same as the ^ = case. 
The last term in Eq. (4.11) is completely negligible after 
{S(/)^) reaches to its maximum value, and the final fiuctu- 
ation can be analyzed in the same way as the minimally 
coupled case. Numerical values of the maximum point 
(,5(/.2)^ = 1.7 X 10-^M|l and M^f) = 1.4 x lO^^ are 
almost the same as those in the minimally coupled case, 
and parametric resonance terminates before the energy of 
the (j)Q field is sufficiently transferred to the inflaton fluc- 
tuation. In this case, the enhanced mode of the momen- 
tum exists also in the narrow range which is estimated 



by Eq. (4.6) 



For C ^ ~ lOi the nonminimal coupling becomes so 
strong that the structure of the resonance changes. Let 
us first consider the dynamics of preheating in the strong 
coupling limit: |^| ^ 1, and later we show the concrete 
cases. 

We assume that the amplitude of the infiaton field is 
not very small, i.e., |^| ^ a^ ;:^ 1 at the initial stage of 
preheating, where a = ^k^0q. When the (j>o field oscil- 
lates far from the minimum of the potential {(f)o ~ 0o) 



we can approximately express the scalar curvature (4.9) 
as 



a \ bt,a 



1 



d<j)o 
dx 



(4.12) 



where (/jq = <i>a / A/pL . Here we did not neglect the term 
l/(6^a) in Eq. ( 4.12| ). This term increases with the pas- 
sage of time because |a| decreases from the order of unity 
as long as (i50^) 
equation of the 

ll 
dx^ 



<.(j)l. We find from Eq. ( [4.1CD that the 
00 field is written by 



1 
6? 



1 



+ 1 /3 « 0. 



(4.13) 



Note that the changing rate of the / field becomes very 
small around (/)o ~ 4>o in the case of the strong coupling. 
This behavior is clearly seen in Fig. 4, which shows the 
evolution of the 4> field in the ^ = —100 case. 

During one oscillation of the inflaton field, |(1 — &^)a\ 
becomes small around the minimum of the potential. In 
this case, the scalar curvature is approximated as 



R w 487rC 



^0 + 4 



dx 



(4.14) 



Then the 00 field satisfies the following equation of mo- 
tion 






i-^s^eui-^ 



dx 



f^O. (4.15) 



In the very vicinity of 0o — 0, Eq. ( 4.15 ) yields 

2' 



dx"^ 



a^ V da; 



f 



0. 



(4.16) 



Because of the existence of the last term, the oscilla- 
tion becomes rapid compared with the minimally coupled 
case (See Fig. 4). During one oscil latio n of t he in flaton 
field, the equa tion changes from ( 4.13 ) to ( 4.15 ), and 
from ( 1.15| ) to ( [4.13D as long as the amplitude of the 0o 
field satisfies the condition a^ ^ 1. With the passage of 
time, however, since the amplitude (po decreases by the 
expansion of the universe, the equation of the 0o field 
approaches to that of the minimally coupled case. 

Next, let us consider the fluctuation of the infiaton 
field. The particle production occurs when the frequency 
( ^.28 ) changes nonadiabatically [|2|. Since the 0o field 
rapidly varies around 0o ~ for the strong coupling 
1^1 ^ 1, the fluctuation grows nonadiabatically wh en 4>q 



is pa ssing through zero. Substituting the relation ( |4.14| ) 
into (4.11), the equation of the fluctuation around (/)o ~ 
is approximately written by |28|] 






<? + 3/2 - A&T^a^i^ 



dx 



Sifik ~ 0. 



(4.17) 



In the vicinity oi 4>q ^ 0, the iSn^"^ ( ^) term dom- 
inates over other terms. We have numerically checked 
that this term changes most significantly around 0o ~ 
and mainly causes the nonadiabatic increase of the fluc- 
tuation. The structure of resonance is modified by this 
effect, and the range of the enhanced momentum modes 
becomes broader. As was presented by Greene et al. Q, 
the minimally coupled case (i.e. the prefactor of /^ is 3) 
is the least favorable case for an efficient preheating. In 
the present model, however, the effect of the nonminimal 
coupling makes the resonance band wider and we can 
expect the rapid growth of the fluctuation. 

Let us consider the case of ^ = —20. In this case, 
initial values for preheating are (f>Q{ti) = 0.0477MpL 
and \ai\ = 1.145 (See Table I). At the initial stage 
of preheating, s ince the structure of resonance is dom- 
inated by Eq. ( 4.17 ) around (Jjq = where the parti- 
cle production occurs, momentum modes close to k — 
can be enhanced. In Fig. 5 we show the evolution of 

wi = 3/2 (l + ^^\ and LU2 = (C - i) a^R terms in 
Eq. ( 4.11| ). We find that the nonminimally coupled term 
u)2 has larger amplitude than uji at the initial stage. 
Since the periods of both functions are the same, UJ2 
dominates over lui for the nonadiabatic change around 



(j>o — 0. However, this term becomes negligible with the 
decrease of the scalar curvature. The increase of {Scj)"^) 
occurs from x « 4000, because the particle creation rate 
can not surpass the expansion rate of the Universe for 
X < 4000. Since the contribution of uj2 is about 10% 
of 071 at a; « 4000, momentum modes close to fc = 
are marginally enhanced when (Scj)^) begins to increase. 
After X w 6000, this contribution becomes less than 
5% and the structure of resonance approaches to the 
minimally coupled case. Actually we can see in Fig. 6 
that the growth rate of {6(j)'^) becomes a bit small for 
X > 6000. Although the range of enhanced k modes is 
wide at the initial stage of preheating, it becomes nar- 
row gradually as preheating proceeds and only the mo- 
mentum modes which satisfy the relation (46) are en- 
hanced. After X « 6000, the fluctuation increases up to 
(5(t)^)f = 1.3 X 10-^M|l, (at Xf = 1.04 x 10*), which is 
almost the same as the ^ = case. Although the reso- 
nant structure is different at the initial stage, this does 
not affect the final abundance of the fluctuation signifi- 
cantly. One important property is that the final ampli- 
tude of the (j)o field is (j)o{xf) — 8.0 x 10^^, which yields 
{Scfy^}/ « 0.2(/)q. This means that the energy of the ho- 
mogeneous inflaton field is more efficiently transferred to 
the inflaton fluctuation than in the minimally coupled 
case. 

With the increase of |^|, the duration of the initial 
resonant stage caused by the nonminimal coupling con- 
tinues longer. However, as time passes and |q;| decreases, 
the structure of resonance approaches the minimally cou- 
pled case. In the case of 1^ — —50, the initial value 
of the inflaton is (t>o{ti) = 0.0303A/pL, which is much 
smaller than that of the ^ = case. Because of this, the 
growth of the inflaton fluctuation is delayed. At the ini- 
tial stage, low momentum modes can be enhanced, and 
this is much more efficient than the minimally coupled 
case with the same initial value of the inflaton. How- 
ever, for 4000 < a; < 7000, {Scj)^) shows a short plateau. 
This means that the range of enhanced k modes becomes 
narrow after x ~ 4000, and the contribution of low mo- 
mentum modes becomes less important. For x > 7000, 
the higher momenta begin to contribute to the growth of 
{5<f)^) ■ Although we can expect the growth of low momen- 
tum modes even for x > 7000, it is rather inefficient and 
the fluctuation produced at the final stage is mainly dom- 
inated by the higher momenta which are estimated by 
Eq. (||6|). The final variance is {5(jP') f = 9.5 x lO'^Ml^ 
at Xf = 2.4 X 10^. The resonance stops by the back re- 
action effect of created particles and the final variance 
does not change so much compared with the ^ = case. 
However, since (l)o{xf) = 4.1 x 10""', we can estimate 
as {S(j>'^)f — 0.5600 ^t the termination point of the res- 
onance. This means that the energy transfer from the 
(j)Q field to 6(l)k particles becomes more efficient with the 
increase of |^|. 

We show in Fig. 7 the evolution of the fluctuation for 
£_ = —70 and —80 cases. The structure of the resonance 



is qualitatively same as the ^ = —50 case. Since the 
^ effect to the fluctuation continues longer with the in- 
crease of 1^1, the fluctuation reaches the first plateau at 
X = 7000 and 8000 for ^ = -70 and -80 cases, respec- 
tively. After a short plateau, the secondary resonant 
stage turns on and ((5(/>^) begins to increase. However, 
at this stage, the growth rate is not so large as the first 
resonant stage, because the enhanced k modes are re- 
stricted in the narrow range. We have found that the 
final variance slowly decreases with the decrease of ^ for 
-80 <£, < - 50. For example, when f = -50, -70, and 
-80, the final variances are {S(j)'^)f/M^^ = 9.5 x 10"**, 
7.8 X 10"^ and 4.7 x 10~^, respectively. As for the ratio 



of 



)f 



to (/)q at the maximum point of the fluctua- 



tion, {S(j)^)f/(t)l = 0.56, 1.15, 1.16 for C = -50, -70, -80, 
respectively. Numerical calculations confirm that reso- 
nance stops when the fluctuation grows up to (50^)/ ~ (fiQ 
for ^ < — 60. Although (f>o{xf) becomes smaller as ^ de- 
creases, the energy of the oscillating inflaton field is trans- 
ferred efficiently to the inflaton particles for ^ < — 60 
as compared with the minimally coupled case. This is 
mainly due to the fact that nonadiabatic creation of the 
inflaton particles efficiently occurs around the minimum 
of the potential. 

For ^ < — 100, the structure of resonance is dominated 
by the ^ effect. We show the evolution of the inflaton 
fluctuation in Fig. 8 for the ^ = -100, -200 and -1000 
cases. We can not separate the resonant stage into two 
parts, but find that low momentum modes are enhanced 
continuously until {6(fP) reaches to its maximum value. 
Because of the nonexistence of the secondary inefficient 
resonant stage, the resonance terminates rather early and 
(l)o{xf) takes larger value than that in the — 80 < ^ < —50 

case. For example, when ^ = —200, 0o ~ 1-0 x 10^^, 
which is much larger than the value (po w 2.0 x 10^** in 
the ^ = —80 case. This increase of the final value of (I>q 
results in the increase of the final fluctuation, since reso- 
nance ends by the back reaction effect when {Scj)"^) grows 
up to the comparable value with (J)q. For example, the 



final values are {S(f>^}f = 2.35 x IO'^M^l for ^ = -100; 
{S(t,'^)j = 1.86 X lO-^A/^L for ^ = -200. For ^ < - 200, 
however, since the final value of (/>o becomes smaller with 
the decrease of ^ again due to the decrease of the ini- 
tial value of (j)o, the final variance also decreases. For 



example. 



1.03 X 10-^Af^L for ^ 



-500, and 



PL for ^ = —1000. In the case of 



\u^2) ^ 5.27 X lO-^Ml 
^ < — 200 the nonadiabatic change of the fluctuation in 
the vicinity of (/)o = is clea rly s een in Fig. 8. This is be- 
cause the ^Ra^ term in Eq. ( 2.28| ) rapidly changes around 
(/>o = especially when |^| is very large as |^| > 200. 
However, in this case also, the growth of the fluctuation 
is finally suppressed by the back reaction effect of created 
particles. 

We have found that the final variance of the nonmin- 
imally coupled case does not change significantly com- 
pared with the minimally coupled case as is found in 



Fig. 9, where the maximum value of {S(fP) is {S(j)^)f = 
2 X 10~^MpL at ^ ~ —200. However, the dispersion of the 
final fluctuation becomes comparable to the amplitude of 
the (f)a field for ^ < - 60 (See Fig. 10). This imphes that 
although the initial value of inflaton decreases as increase 
of 1^1, this is compensated by the efficient resonance with 
the appearance of new resonant modes close to fc = 0, We 
depict in Fig. 11 the distribution of produced momentum 
modes for the ^ = 0, —1, —50, —100 cases. With the de- 
crease of ^, we find that momentum modes close to fc = 
are mainly enhanced. When £, < — 100, the final fluc- 
tuation is dominated by low momentum modes, and this 
is different from the minimally coupled case. As we have 
seen, preheating stage generically exists even if the non- 
minimal coupling is taken into account. 

Finally, we should mention the case when A is changed. 
For largely negative ^, larger values of A is allowed to fit 
t he density perturbation because the effective potential 
( ^.7[ ) be com es flat with the increase of |^|. As is found 
by Eq. (4.4), the frequency of the 0o field becomes large 
with the increase of A, which results in the fast growth 
of the fluctuation. However, A is scaled out in the nor- 
malized equations of the 00 and S4>k fields. A appears 
only in the regularization of the initial amplitude of the 
fluctuation. Since the growth of the fluctuation stops 
by the back reaction effect of created particles, the final 
variance is almost the same as the case of the fixed A, as 
was suggested in the minimally coupled case in Ref . [Q . 



V. CONCLUDING REMARKS AND 
DISCUSSIONS 

In this paper, we have investigated preheating of the 
nonminimally coupled inflaton field in the chaotic infla- 
tion model. As for the massive inflaton case, since ^ is 
constrained to be small as |^| < 10"'^ to lead a sufficient 
inflation, the evolution of the system is almost the same 
as the minimally coupled case. In the self-coupling in- 
flaton case, the constraint of ^ is absent for negative £,. 
Hence we considered mainly this case and investigated 
how the ordinary picture of preheating is modified by 
taking into account the ^ effect. 

In the minimally coupled case, the equation of the in- 
flaton fluctuation is reduced to the Lame equation at the 
linear stage of preheating. Although this model is not 
so efficient for the development of the resonance com- 
pared with the massive inflaton plus another scalar field 
X coupled to inflaton, the fluctuation {6<j)'^) increases by 
parametric resonance and reaches to its maximum value 
((5(/)2) = 1.5 X 10~^MpL. In reheating phase, the time- 
averaged scalar curvature vanishes in this model, and the 
investigation of preheating is reduced to the theory of 
Minkowski spacetime. Since the instability band is nar- 
row as 3/^/2 < p^ < V?>p , the resonance is sensitive to 
the back reaction effect of created particles. The growth 
of the fluctuation ends at {Scfy^)/ « O.O50§, which means 



that the dispersion of the fluctuation is only 20% of the 
amplitude of the homogeneous inflaton field. 

Taking into account the nonminimal coupling with the 
scalar curvature, the dynamics of preheating can be al- 
tered. We found that the i^ effect works in two ways. 
First, since the initial value of inflaton for preheating 
becomes lowered as |^| increases, and the frequency of 
the oscillation of the inflaton field is reduced. Hence the 
growth of the inflaton fluctuation is delayed by this ef- 
fect for fixed A. Second, since frequencies of the both 
homogeneous and fluctuational parts change due to the 
nonminimal coupling, the structure of resonance becomes 
different from the minimally coupled case. In the case of 
—O.I < ^ < 0, the evolution of the system is almost the 
same as the minimally coupled case. In —10 ^ C ^ ~ 0-1 
case, the decrease of the initial value of inflaton mainly 
delays the evolution of the system. Although the £, ef- 
fect changes both frequencies of the (po and S(f>k fields 
at the initial stage of preheating, this works only for a 
short period. When the particle creation rate surpasses 
the expansion rate of the Universe and {S(p'^) begins to 
increase, the £, effect on the (f>o and 5(pk fields can be 
negligible. Hence, after that, the system of the evolution 
is not affected by the ^ effect. The structure of reso- 
nance seems almost the same as the minimally coupled 
case, and produced momentum modes are restricted in a 
narrow range when (Scj)^) reaches to its maximum value. 
The final variance does not change significantly compared 
with the minimally coupled case. 

As £, decreases further {£, < — 10), the initial stage of 
preheating is dominated by the effect of the nonminimal 
coupling. The structure of resonance becomes different 
from that of the ^ = case, and the wider range of mo- 
menta is allowed. Particle creation occurs nonadiabati- 
cally in the vicinity of 4>o = 0. Although the minimally 
coupled case is the least favorable one for the develop- 
ment of the resonance because of the limited resonant 
band, such limitation can be removed by taking into ac- 
count the nonminimal coupling. When — 80<^< — 10, 
the stage of parametric resonance is separated into two 
stages. At the first stage, the resonance caused by the 
nonminimal coupling dominates, and momentum modes 
close to A: = are efficiently enhanced. With the passage 
of time, the S, effect onto equations of the ^o a-nd S(j)k 
fields becomes less important and the system enters the 
secondary resonant stage. The second stage is the ordi- 
nary resonance stage of the minimally coupled case, and 
the growth rate of the fiuctuation is not so large as that 
of the first stage. It is true that low momentum modes 
are enhanced, but the final fluct uati on is dominated by 
higher modes described by Eq. (|4.6|) . Although the fi- 
nal variance slowly decreases with the decrease of £, for 
—80 < ^ < ^10, the ratio of the final fiuctuation {5(j)^) f 
to the square of the amplitude of the (/)o field gets larger. 
Especially for S, < — 60, we find {64>'^)f w 0§, which im- 
plies that the transfer of the energy from the 4>o field to 
the fluctuation is much more efficient than the minimally 



coupled case. 

When C < — 100, the growth of fluctuation ends before 
the secondary resonant stage appears, and parametric 
resonance is quite efficient. For —200 ^ £, ^ — 100, 
the final variance increases with the decrease of ^ and 

{5(1)^) f takes the maximum value {S<j)'^)f = 2 x IO^^A/|t 

a2\ 



at ^ « —200. For ^ < — 200, however, {Scf)"^)/ decreases 
because the initial value of infiaton becomes smaller with 
the decrease of £,, which yields the smaller amplitude at 
the termination point of resonance. We also found the 
relation (<50^)/ « <^o for ^ < — 100, which means that 
resonance mainly terminates by the back reaction effect 
when the dispersion of the fluctuation grows up to the 
comparable value of the amplitude of the ^o field. 

Finally we have to comment on some points. In this 
paper, we have made heavy use of the Hartree approxi- 
mation. In taking into account the back reaction effect 
of created particles, there are other approaches. One 
of them is to add the stochastic noise term originated 
by quantum fluctuation to the inflaton equation. This 
is based on the closed time path formalism in the non- 
equilibrium quantum field theory p3|, and several au- 
thors made use of this method in the context of preheat- 
ing pOJ . It is expected that the growth rate of the fiuc- 
tuation would be changed by taking into account this ef- 
fect. Other approach is the lattice simulation. Although 
the final fiuctuation of our numerical calculation in the 
minimally coupled case coincides well with that of the 
full nonlinear calculation in Ref. [|3|, the final variance 
may change in the nonminimally coupled case performed 
by the lattice simulation. Since the lattice simulation 
includes the rescattering effect which is absent in the 
Hartree approximation, the spectrum of enhanced mo- 
mentum modes would also change by rescattering of cre- 
ated particles. Although we do not present details of 
the difference in this paper, we should consider preheat- 
ing by various approaches of the back reaction effect and 
compare them with the fully nonlinear calculation. 

Recently, several authors considered the metric pertur- 
bation in reheating phase [g^. As for the self-coupling 
inflaton model, it was pointed out in Ref. ||3^ that the 
metric perturbation is resonantly amplified even if we 
do not introduce another scalar field coupled to infia- 
ton. In the present model of the nonminimal coupling, 
since the super Horizon modes can be enhanced by the 
strong nonminimal coupling, this may alter the spectrum 
of the density perturbation. It will be interesting to in- 
vestigate preheating including the metric perturbation, 
because modification of the spectrum of the density per- 
turbation would result in some important consequences 
such as the production of primordial black holes. These 
issues are under consideration. 



10 



ACKOWLED CEMENTS 

S. T. would like to thank B. A. Bassett and A. 
Taruya for useful discussions. T. T. is thankful for fi- 
nancial support from the JSPS. This work was sup- 
ported partially by a Grant-in- Aid for Scientific Research 
Fund of the Ministry of Education, Science and Cul- 
ture (No. 09410217 and Specially Promoted Research 
No. 08102010), by a JSPS Grant-in- Aid (No. 094162), 
and by the Waseda University Grant for Special Research 
Projects. 



[1] J. Traschen and R. H. Brandenberger, Phys. Rev. D 42, 
2491 (1990); Y. Shatanov, J. Trashen, and R. H. Bran- 
denberger, Phys. Rev. D 51, 5438 (1995). 

[2] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. 
Lett. 73, 3195 (1994). 

[3] A. D. Dolgov and A. D. Linde, Phys. Lett. B116, 329 
(1982); L. F. Abbott, E.Fahri, and M. Wise, Phys. Lett. 
B117, 29 (1982). 

[4] E. W. Kolb, A. D. Linde, and A. Riotto, Phys. Rev. 
Lett. 77, 3716 (1996); E. W. Kolb, A. Riotto, and L L 
Tkachev, Phys. Lett. B423, 348 (1998). 

[5] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. 
Lett. 76, 1011 (1996); L L Tkachev, Phys. Lett. B376, 
35 (1996); S. Kasuya and M. Kawasaki, Phys. Rev. D 56, 
7597 (1997). 

[6] S. Khlebnikov, L. Kofman, A. Linde, and L I. Tkachev, 
Phys. Rev. Lett. 81, 2012 (1998); L L Tkachev, S. Khleb- 
nikov, L. Kofman, A. Linde, Phys. Lett. B440, 262 
(1998); S. Kasuya and M. Kawasaki, Phys. Rev. D 58, 
083516 (1998). 

[7] S. Khlebnikov and L L Tkachev, Phys. Rev. D 56, 653 
(1997). 

[8] B. A. Bassett, Phys. Rev. D 56, 3429 (1997). 

[9] S. Khlebnikov and L L Tkachev, Phys. Lett. B390, 80 
(1997). 
[10] S. Khlebnikov and I. L Tkachev, Phys. Rev. Lett. 79, 

1607 (1997). 
[11] T. Prokopec and T. G. Roos, Phys. Rev. D 55, 3768 
(1997); B. R. Greene, T. Prokopec, and T. G. Roos, Phys. 
Rev. D 58, 6484 (1997). 
[12] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. 

D 73, 3258 (1997). 
[13] S. Khlebnikov and I. L Tkachev, Phys. Rev. Lett. 77, 

219 (1996). 
[14] D. Boyanovsky, H. J. de Vega, R. Holman, D. S. Lee, 
and A. Singh, Phys. Rev. D 51, 4419 (1995); D. Boy- 
anovsky, M. D' Attanasio, H. J. de Vega, R. Holman, D. 
S. Lee, and A. Singh, Phys. Rev. D 52, 6805 (1995); D. 
Boyanovsky, H. J. de Vega, R. Holman, D. S. Lee, and A. 
Singh, J. F. J. Salgado, Phys. Rev. D 54, 7570 (1996); D. 
Boyanovsky, D. Cormier, H. J. de Vega, R. Holman, A. 
Singh, and M. Srednicki, Phys. Rev. D 56, 1939 (1997); 



D. Boyanovsky, R. Holman, S. Prem Kumar, Phys. Rev. 

D 56, 1958 (1997). 
[15] J. Baacke, K. Heitmann, and C. Patzold, Phys. Rev. D 

55, 2320 (1 997); ibid D 56, 6556 (1997); J. Baacke and 

C. Patzold, [hep-ph/9906417 . 
[16] D. L Kaiser, Phys. Rev. D 53, 1776 (1995); ibid D 56, 

706 (1997). 
[17] D. T. Son, Phys. Rev. D 54, 3745 (1996). 
[18] P. B. Greene, L. Kofman, A. Linde, and A. A. Starobin- 
sky, Phys. Rev. D 56, 6175 (1997). 
[19] B. A. Bassett and S. Liberati, Phys. Rev. D 58, 021302 

(1998). 
[20] S. Tsujikawa, K. Maeda, and T. Torn, Phys. Rev. D 60, 

063515 (1999). 



[21] S. Tsujikawa, K. Maeda, and T. Torii, tiep-ph/9906501 , 
to appear in Phys. Rev. D. 

[22] T. Futamase and K. Maeda, Phys, Rev. D 39, 399 (1989). 

[23] R. Fakh and W. G. Unruh, Phys, Rev. D 41, 1783 (1990). 

[24] N. Makino and M. Sasaki, Prog. Theor. Phys. 86, 103 
(1991). 

[25] E. Komatsu and T. Futamase, Phys, Rev. D 58, 023004 
(1998). 

[26] In other papers, the value of A = 10^^'^ is often adopted. 
However, we have confirmed that the difference for these 
two cases of A is negligibly small and hardly alters our 
results and discussions. 

[27] H. Kodama and T. Hamazaki, Prog. Theor. Phys. 96, 949 
(1996); Y. Nambu and A. Taruya, Prog. Theor. Phys. 97, 
83 (1997); B. A. Bassett, D. L Kaiser, and R. Maartens, 
Phys. Lett. B455, 84 (1999); B. A. Bassett, F. Ta m- 



burini, D. L Kaiser, and R. Maartens, liep-ph/990131£ , 
to appear in Nucl. Phys. B; F. Finelli and R. Branden- 
berger, Phys. Rev. Lett. 82 1362 (1999); M. Parry and 
R. Easther, Phys. Re v. D 59 061301 (1999); B. A. Bas- 
sett and F. Viniegra, ^iep-ph/9909353| ; B. A. Bassett, C . 
Gordon, R. Maartens, and D. L Kaiser |hep-ph/9909482 . 



[28] We have also obtained the approximate equation of 
the fluctuation around <jf>o — (j>o- This is expressed as 
-j-^Sipk + (p^ + 2/^j Sifk ~ 0, and the structure of reso- 
nance is difi^erent from that of the minimally coupled case. 
In this case, however, since the (j)o field moves very slowly, 
the nonadiabatic increase of the fiuctuation hardly oc- 
curs. 

[29] E. Calzetta, A. Campos, and E. Verdaguer, Phys. Rev. 
D 56 2163 (1997); E. Calzetta and B. L. Hu, Phys. Rev. 
D 52 6770 (1995); A. Campos and B. L. Hu, Phys. Rev. 
D 58, 125021 (1998). 

[30] S. A. Ramsey, B. L. Hu, and A. M. Stylianopoulos, Phys. 
Rev. D 57, 6003 (1998); B. A. Bassett and F. Tam- 
burini, Phys. Rev. Lett. 81, 2630 (1998); B. A. Bassett, 
Phys. Rev. D 58, 021303 (1998); V. Zanchin, A. Maia 
Jr., W. Craig, and R. Brandenberger, Phys. Rev. D 57 
4651 (1998); V. Zanchin, A. Maia Jr., W. Craig, and R. 
Brandenberger, Phys. Rev D 60, 023505 ( 1999). 

[31] M. Parry and R. Easther, |hep-ph/990355C . 



11 



TABLE I. The final value of |a| and (f>o when the infla- 
tionary period ends in the A<^'*/4 model for negative ^. We 
find that lap] takes the almost constant value |q;f| ~ 1 for 
5 < —1. On the other hand, (jioitp) becomes smaller as |5| 
increases. 
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Figure Captions 

FIG. 1: 

The potential in the equivalent system to the non- 
niinimaUy coupled inflaton. The dotted and solid 
curves denote the massive inflaton and the massless 
inflaton case with ^ = —0.1 respectively. The po- 
tential of the massive inflaton has a local maximum 
and inflation is difficult to realize unless |^| < 10^'^. 
However, the potential of the massless inflaton has a 
flat plateau to lead a sufficient inflation for negative 

i- 

FIG. 2: 

The evolution of {5(p) as a function of a; in ^ = 
and ^ = — 1 cases. In both cases, the fluctuation 
grows by parametric resonance. When ^ = — 1, since 
the initial value of inflaton becomes smaller com- 
pared with the ^ = case, this makes the growth of 
the fluctuation delayed. However, the final variance 
{Sep) ~ 10~^ is almost the same as the minimally 
coupled case. 

FIG. 3: 

The evolution of the homogeneous inflaton field as a 
function of a; in ^ = and ^ = — 1 cases. Since the 
initial value of inflaton for the ^ = — 1 case is smaller 
than that of the ^ = case, the time interval of the 
oscillation becomes longer. 

FIG. 4: 

The evolution of the (po field as a function of x in 
the case of ^ — —100. The 0o field slowly changes 
around the local maximum point of the potential, but 
nonadiabatically changes in the vicinity of <f>Q w 0. 

FIG. 5: 



The evolution of ui = 3p 1 



p 



and 



UJ2 



(C ~ ^) '^^R ^s a function of x in the case of ^ == 
—20. The effect by the W2 term is significant at the 
initial stage of preheating, but it gradually becomes 
negligible with the passage of time. 

FIG. 6: 

The evolution of {S(fP) as a function of a; in ^ = -20 
and £, = —50 cases. At the initial stage of preheating, 
the structure of resonance is modified and the range 
of enhanced momenta becomes wider. However, this 
stage ends before (ScjP) sufficiently increases and the 
second resonance stage sets in. In both cases of ^ = 
-20 and £_ = —50, the growth rate of the fluctuation 
becomes smaller after the first resonance stage ends. 

FIG. 7: 
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The evolution of {5(p^) as a function of x in ^ = — 70 
and ^ = —80 cases. We can separate the resonance 
into two stages. The first stage is dominated by the 
effect of the nonminimal couphng, which connects to 
the second stage of ordinary minimal coupling. With 
the decrease of ^, the initial stage becomes longer. 

FIG. 8: 

The evolution of ((50^) as a function of a; in ^ = —100, 
^ = —200 and ^ = —1000 cases. In these cases, the 
fluctuation reaches to its maximum value only by the 
first efficient resonance stage by the effect of the non- 
minimal coupling. Therefore, the final fluctuation is 
dominated by the momenta close to fc = 0. 

FIG. 9: 

The final value of {5(IP) as a function of |^| for the 
negative ^. In the case of — 80 < ^ < 0, {S4>'^)f slowly 
decreases with the increase of |^|. For ^ < — 100, 
the resonance becomes efficient by the ^ effect, and 
{Scjr)/ takes the maximum value {S4>'^)f « 2 x 10~^ 
at ^ w —200. However, {5(p) f slowly decreases for 
^ < — 200 because the growth of the fluctuation 
is suppressed by the back reaction effect of created 
particles. 

FIG. 10: 

The ratio of the final variance (50^)/ to the square 
of the amplitude of the ijio field. This ratio increases 
with the decrease of ^ for — 50 < f < 0, and reaches 

the plateau {S(f>'^)f/4>1 w 1 for ^ < - 60. 

FIG. 11: The spectrum of enhanced momentum modes in 

C = 0, C = -1, C = -50, and ^ = -100 cases. 
Although the resonance band is narrow in the min- 
imally coupled case, the low momentum modes be- 
gin to be enhanced with the decrease of ^. When 
^ < — 100, the final fluctuation is dominated by low 
momentum modes. 
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